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Abstract. Warm Dark Matter (WDM) models have recently been resurrected to 
resolve apparent conflicts of Cold Dark Matter (DM) models with observations. En- 
dowing the DM particles with non-negligible velocities causes free-streaming, which 
suppresses the primordial power spectrum on small scales. The choice of a root-mean- 
square velocity dispersion Urms,o ^ 0.05 km/s at redshift z = Q (corresponding to a 
particle mass mx ~ 1 keV if the WDM particles are fermions decoupling while rela- 
tivistic) helps alleviate most, but probably not all, of the small-scale problems faced 
by CDM. An important side-effect of the particle velocities is the severe decrease in 
the number of collapsed halos at high redshift. This is caused both by the loss of 
small-scale power, and by the delay in the collapse of the smallest individual halos 
(with masses near the effective Jeans mass of the DM). The presence of early halos 
is required in order (1) to host either early quasars or galaxies that can reionize the 
universe by redshift z = 5.8, and (2) to allow the growth of the supermassive black 
hole believed to power the recently discovered quasar SDSS 1044-1215 at this redshift. 
We quantify these constraints using a modified Press-Schechter formalism, and find 
Wrms,o ^ 0.04 km/s (or mx ^ 1 keV). If future observations uncover massive black 
holes at z > 10, or reveal that reionization occurred at 2 ^ 10, this could conclusively 
rule out WDM models as the solution to the small-scale crisis of the CDM paradigm. 



INTRODUCTION 



The currently favored model of hierarchical galaxy formation in a universe dom- 
inated by cold dark matter (CDM) has been very successful in matching observa- 
tions of the density distribution on large scales. However, recently some small-scale 
shortcomings of this model have appeared, as summarized elsewhere in these pro- 
ceedings (see also [1] for a recent review). Although the observational significance 
of the discrepancies is still disputed, and astrophysical solutions involving feedback 
may still be possible, the accumulating tension with observations has focused atten- 
tion on solutions involving the particle properties of dark matter. Proposals have 
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included self-interacting dark matter [2], adding a repulsive interaction to gravity 
[3,4], the quantum-mechanical wave properties of ultra-light dark matter particles 

[5], and a resurrection of warm dark matter (WDM) models [6-8]. By design, a 
common feature of models that attempt to solve the apparent small-scale prob- 
lems of CDM is the reduction of fluctuation power on small scales. In the CDM 
paradigm, structure formation proceeds bottom-up: the smallest objects collapse 
first, and they subsequently merge together to form larger objects. It then follows 
that the loss of small-scale power modifies structure formation most severely at 
the highest redshifts; in particular, the number of self-gravitating objects at high 
redshift is reduced. 

A strong reduction in the abundance of high-redshift objects could be in conflict 

with observations. First, the lack of a Gunn-Peterson trough in the spectrum of 
the bright quasar SDSS 1044-1215 at redshift z = 5.8 [9] implies that the hydrogen 
in the intergalactic medium (IGM) was highly ionized prior to this redshift. The 
most natural explanation for reionization is photo-ionizing radiation produced by 
an early generation of stars or quasars [10]. The sources of reionization reside in 
halos that have masses in the range corresponding to dwarf galaxies — the mass 
scale on which power needs to be reduced relative to CDM models. Second, if SDSS 
1044-1215 is unlensed and radiating at or below the Eddington limit, its unusual 
intrinsic brightness implies that it is powered by an exceptionally massive black 
hole (BH). The growth of this BH, out of a stellar remnant seed, requires [11] that 
a host halo be present at a sufficiently high redshift {z ^ 5.8). 

In this contribution, we briefly review the status of WDM models in resolving 
the problems of CDM, and then examine new constraints that arise on WDM 
models from the high redshift universe. We focus on WDM models, although similar 
constraints would apply to other modifications of the CDM paradigm that reduce 
the small-scale power. The cosmological parameters we adopt, based on present 
large scale structure data [12], are (fio, ^a, f^b, h, as, n) = (0.3, 0.7, 0.045, 0.7, 0.9, 1). 
The details of the study described here can be found in [13]. 

WARM DARK MATTER MODELS 

The WDM is assumed to be composed of particles of about ~ 1 keV mass (com- 
pared to ~ 1 GeV in CDM, or ~ 10 eV in Hot Dark Matter models). The thermal 
velocities of the particles cause free streaming out of overdense regions, smoothing 
out small-scale fluctuations, leading to a small-scale cutoff in the linear power spec- 
trum. In addition, the thermal velocities act similarly to pressure, and inhibit the 
growth of low-mass perturbations. One example of WDM is fermionic particles 
that decouple in the early universe while relativistic and in thermal equilibrium 
[14]. To produce a given contribution Vlx to the cosmological critical density, the 
required particle mass mx is then determined by mxnx oc ^xh^, where the present 
number density nx of WDM particles follows from their chosen r.m.s. velocities. 
This yields a relation between particle mass and r.m.s. velocity dispersion [8], 



where gx is the effective number of degrees of freedom of WDM. The comoving 
cutoff scale i?c, where free-streaming reduces the power spectrum to half of its 
value in CDM, is given by 

/O h2 \ -0-14 / \ 0.86 

= 0.226(^1 ( ^-7'° , ) Mpc. (2) 
\^ 0.15 I VO-05 km/sy ^ ^ ^ 

The length scale Rc corresponds to a characteristic mass scale Mc, 

M, = 1.74 X 10« f f Mo . (3) 

The predictions of WDM models can be expected to differ from CDM on scales 
below Rc or Mf.. In addition, it is useful to define an "effective Jeans mass" for 
WDM, at which the pressure corresponding to the r.m.s. particle velocities balances 
gravity, 

M, = 3.06 X 10» f^V' (i±ii)"' Me 

Vi.5/ V 0-15 y Vikev; Vsooo; ^ 

= 4.58 X 10» (^Y" ( V, fi±ii)'" . (4) 
V 0.15 I km/s y V 3000 J ^ ^ ' 

We have verified these scalings using spherically symmetric, one-dimensional col- 
lapse simulations [13]. The growth of perturbations on scales below Mj is slowed 
down by the "pressure" of WDM. Although this effect is irrelevant in most discus- 
sions of WDM models at lower redshifts (due to the smallness of Mj), we find that 
it is important in the context of reionization. The inhibited growth of perturba- 
tions in the linear regime results in a delay in the final virialization epoch of the 
individual, low-mass halos, that first condense out in the universe. 



SUMMARY OF CURRENT OBSERVATIONS 

The current status of WDM models in resolving various problems of CDM is 
summarized in Table 1. Based on analytical arguments, as well as numerical sim- 
ulations, WDM particles with a mass in the range 0.6 keV ^ mx ^ 1-5 keV 

can strongly suppress the number of satellite halos of the Milky Way, and produce 
low halo concentration parameters for both dwarfs, and Milky Way-sized galaxies 



TABLE 1. Scorocard of WDM models. 



CDM Problem'' WDM mass Status References 



Milky Way satellites 


0.6- 


-1.5 keV 


V 


[6,8,16] 


Halo concentration 


0.6- 


-1.5 keV 


V 


[6,8,17] 


Inner density profile 


<0.6 keV ? 




[6,17] 


No dwarfs in voids 


<1A 


) keV 


? 


[8,19] 


Angular momentum 


0.5- 


-0.8 kcV 


X? 


[7,20] 



^ Note that WDM particles with nix ^ 750 eV would 
contradict the power spectrum of the Lya forest [15]. A 
similar limit follows from the maximum observed phase 
space density in cores of dwarf spheroidals [17,18]. 



[6,8]. It remains unclear whether WDM can fully resolve some other, related prob- 
lems. One of the major successes of the CDM paradigm is the interpretation of 
the statistics of the Lya forest. An analysis of the Lya forest power spectrum at 
redshift z ~ 3 finds that WDM models with mx < 0.75 kcV would spoil this suc- 
cess [15] . A similar limit from the maximum observed phase space density in dwarf 
spheroidal galaxies. To solve some of the problems with CDM, lower masses might 
be needed. In existing simulations for rux > 0.6 keV, the inner slope of the density 
profiles of normal disk galaxies do not appear to fiatten sufficiently to produce the 
observed, nearly constant cores, although these simulations do not yet probe the 
relevant innermost regions reliably [6]. A similar situation arises with two other 
CDM problems. The solution of the long-standing angular-momentum problem 
in galaxy formation might be solved if "subclumps" in the protogalaxy's halo were 
eliminated (CDM models produce too small disks, thought to be attributable to an- 
gular momentum transfer from the infalling gas to the DM subclumps). Although 
WDM does suppress halo substructure, the required WDM mass is estimated to 
be 0.5 — 0.8 keV [7]. Furthermore, a simulation that artificially turns off angular 
momentum transfer between gas and DM halo subclumps reveals that the specific 
angular momentum distribution still contradicts observations [20] . Finally, the ab- 
sence of dwarf galaxies in voids, a difficulty of CDM [19], is found to be eased in 
a large-scale simulation of an mx = 1.5 keV WDM model [8], but a smaller mx 
appears necessary to match the observations in detail. 

CONSTRAINTS FROM REIONIZATION 

As emphasized in section § 1, WDM models reduce the abundance of high- 
redshift halos, which could lead to confiicts with observations. Note that in general, 
the differences between WDM and CDM models are amplified at high redshifts. In 
a study based on a modified version of the extended Press-Schechter theory, we 
have quantified constraint on the velocities and masses of WDM particles from 
reionization, and from the presence of a supermassive BH in the bright z — 5.8 
quasar SDSS 1044-1215. The main features of our models are as follows: 



1. The effect of free streaming on the power spectrum was included using a fit to 
the numerical transfer function, obtained from a Boltzmann code. 

2. The effect of WDM velocities on the growth of perturbations was included 
using spherically symmetric, one-dimensional hydrodynamical simulations. 

3. The WDM halo mass function was computed using Monte Carlo-generated 
halo "merger trees". A moving barrier {6c > 1.69) was adopted to generalize 
(an improved) Press-Schechter theory to the case of mass-dependent collapse 
times. Our semi-analytic halo mass function was demonstrated to agree with 
numerical simulations at low redshift over the halo masses of interest [8]. 

4. We assume that a fraction of baryons turns into stars in halos with virial 
temperatures Tvir > lO^K (necessary for efficient cooling [21]). We assume 
a Scalo (1998) stellar IMF (producing ^ 4000 ionizing photons per baryon), 
and that a fraction /esc of the ionizing photons escape into the IGM. We 
parameterize our models by the product = /*/esc(= 0.01 in our standard 
model), consistent with the ionizing background inferred from the proximity 
effect at redshift z^3. [22] 
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FIGURE 1. The filling factor of ionized hydrogen, Fhii, as a function of redshift z in our 
standard model (C = 10, = 0.01). The uppermost curve corresponds to ACDM, and the lowest 
curve to a WDM particle mass of mx = 1-25 keV. The middle pair of curves shows separate 
contributions to the delay in reionization from the suppression of the power spectrum (dotted 
curve) and the "effective pressure" of WDM (long dashed curve). 



5. The mean clumping of the ionized gas in the IGM is parameterized by a single 
constant C = (ti^) /n'^{= 10 in our standard model). 

Our main results on the constraints from reionization are shown in Figures 1 and 
2. Figure 1 shows the filling factor of ionized hydrogen in our standard model, with 
nix = 1-25 keV, = 0.01, and C = 10. In this model, reionization occurs at ^ = 
5.8. The contributions to the decrease of the filling factor from the suppression of 
the power spectrum and the effective pressure of WDM are comparable (dotted and 
long-dashed curves, respectively). In Figure 2, we demonstrate how the reionization 
redshift depends on nix, and on our two model parameters, e,, and C. 

CONCLUSIONS 

If high-redshift galaxies produce ionizing photons with an efficiency similar to 
their z = 3 counterparts (e* ~ 0.01), reionization by redshift z = 5.8 places a 
limit of mx ^ 1-25 keV (frms,o ^ 0.03 km/s) on the mass of the WDM particles. 
This limit is somewhat stronger than the limit inferred from the statistics of the 
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FIGURE 2. Reionization redshift as a function of WDM particle mass. Our standard model 
for the star-formation efficiency, escape fraction, and stellar IMF, is shown by the thick solid 
line. The three solid curves show models with the same clumping factor (C = 10) but different 
ionization efficiencies (e, = 0.01,0.1,1, bottom to top). The three dotted curves correspond to 
the same models except with C = 1 ; the three dashed curves correspond to the same models with 
C = 30. In our standard model, the requirement 2;icion > 5.8 yields the constraint mx > 1-25 keV. 



Lya forest (which yields nix ^ 0.75 keV; [15]), although our limit may weaken to 
nix ^ 0.75 keV (t'rms.o = 0.060 km/s) given the uncertainty in current measure- 
ments of the stellar contribution to the ionizing intensity at 2; = 3 [22] (which we 
can use to normalize our models). We also find that the existence of a 4 x 10^ Mq 
supermassive black hole ai z = 5.8, powering the quasar SDSS 1044-1215 (assuming 
it is unlensed and radiating at or below the Eddington limit), yields the somewhat 
weaker, but independent constraint mx ^ 0.5 kcV (or t'lms.o ^ 0.10 km/s), if this 
BH acquired most of its mass accreting at the Eddington rate. Finally, we also find 
that WDM models with mx ^ 1 keV (vrms.o 0.04 km/s) produce a low- luminosity 
cutoff in the high- 2; galaxy luminosity function which is detectable with the Next 
Generation Space Telescope. Such an observation would directly break the degen- 
eracy in the reionization redshift between low ionizing-photon production efficiency 
and small WDM particle mass. The constraints derived here will tighten consider- 
ably as observations probe still higher redshifts, offering increasingly stringent tests 
of models with diminished small-scale power, exemphfied by WDM. 

ZH acknowledges support from a Hubble Fellowship, and RB acknowledges sup- 
port from CITA and from Institute Funds (IAS, Princeton). 
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